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INTRODUCTION  AND  SUMMARY 

The  work  reported  here, — funded  under  Contract  N00014-85-C-0765, 
considers  the  computation  of  the  microwave  brightness  temperatures  of  various 
types  of  sea  ice.  The  methodology  is  an  application  of  theoretical  and 
experimental  results  obtained  at  Aerojet  ElectroSystems  Company  under  pre¬ 
vious  contracts  with  the  Office  of  Naval  Research.  In  these  studies,  the 
bilocal  approximation  of  strong  fluctuation  theory  was  used  to  develop 
equations  for  the  effective  dielectric  constant  of  both  isotropic  and  aniso¬ 
tropic  random  media  as  well  as  equations  for  calculating  the  brightness  tem¬ 
perature  of  isothermal  isotropic  media.  The  earlier  work  also  examined  the 
present  state  of  knowledge  concerning  the  dielectric  properties  of  the  con¬ 
stituents  of  sea  ice.  It  was  found  that  there  were  significant  uncertainties 
relating  to  the  loss  tangent  of  pure  ice  as  well  as  the  dielectric  properties 
of  brine  in  sea  ice.  The  latter  problem  was  overcome  by  a  systematic  meas¬ 
urement  of  the  dielectric  properties  of  brine  and  the  derivation  of  Debye 
parameters  which  allow  the  dielectric  constant  of  brine  to  be  calculated 
throughout  the  microwave  portion  of  the  electromagnetic  spectrum.  The  uncer¬ 
tainties  with  respect  to  the  loss  tangent  of  pure  ice  remain  so  that  a  tenta¬ 
tive  model  based  on  available  experimental  data  must  be  used. 

The  specific  tasks  performed  under  the  present  contract  consisted 
of  (1)  the  development  of  a  computer  program  to  calculate  the  brightness 
temperature  of  an  isothermal  layered  structure,  each  of  whose  layers  could 
consist  of  an  anisotropic  random  medium  and  (2)  the  application  of  the  model 
to  a  variety  of  sea  ice  types  at  different  temperatures  and  the  comparison 
of  calculations  with  some  reported  measurements.  Section  1  of  this  report 

ii 
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discusses  the  basic  equations  required  for  the  computer  program  and  concen¬ 
trates  specifically  on  the  modification  of  previously  reported  equations 
demanded  by  the  introduction  of  elongated  brine  cells  leading  to  a  descrip¬ 
tion  in  terms  of  a  tensor  correlation  function.  Section  2  of  this  report 
discusses  computed  brightness  temperatures  of  various  types  of  sea  ice.  In 
order  to  facilitate  comparison  with  experiments,  data  on  new  ice,  young  ice, 
first  year  ice  and  multiyear  ice  are  treated  separately.  It  is  found  that 
the  main  features  of  the  brightness  temperature  measurements  on  the  first 
three  ice  types  can  be  explained  by  the  model  over  a  wide  range  of  tempera¬ 
tures  and  snow  cover  conditions.  For  multiyear  ice,  however,  discrepancies 
between  the  theory  and  experiments  begin  to  appear  somewhere  in  the  20-30  GHz 
range  and  become  progressively  more  pronounced  at  higher  frequencies.  These 
appear  to  be  the  result  of  the  model  underestimating  scattering  effects  on 
the  brightness  temperature  for  multiyear  ice  and  indicates  the  need  for 
further  theoretical  developments.  Finally,  Section  3  contains  conclusions 


and  recommendations. 
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Section  1 

BRIGHTNESS  TEMPERATURE  OP  A  RANDOM  ANISOTROPIC  MEDIUM 

Equations  for  calculating  the  brightness  temperature  of  an  isother¬ 
mal  layered  medium  have  been  discussed  in  [1]  for  the  case  where  the  effective 
dielectric  properties  of  each  of  the  layers  is  characterized  by  a  scalar  die¬ 
lectric  constant  and  correlation  function.  Since  it  has  been  shown  in  [2] 
that  certain  forms  of  sea  ice  have  anisotropic  dielectric  properties  which 
require  a  description  by  means  of  a  second  rank  dielectric  tensor,  it  is  of 
some  interest  to  see  how  the  basic  equations  may  be  generalized  to  the  aniso¬ 
tropic  situation.  This  problem  is  treated  here. 

1.1  RELATION  OF  BRIGHTNESS  TEMPERATURE  TO  SCATTERING  COEFFICIENTS 

For  an  isothermal  medium  there  is  a  close  relationship  between  the 
brightness  temperature  and  bistatic  scattering  coefficients,  which  are  pro¬ 
portional  to  the  mean  intensity  of  the  scattered  field,  when  the  medium  is 
illuminated  by  a  plane  electromagnetic  wave.  Using  the  relationships  devel¬ 
oped  by  Peake  [3]  and  expressing  the  scattering  coefficients  of  a  random 
medium  bounded  by  flat,  parallel  surfaces  as  the  sum  of  a  specular  term  and  a 
bistatic  term  arising  from  the  random  field  within  the  medium,  it  was  shown 
in  [1]  that  the  emissivity  of  the  medium  in  the  direction  -k0  (see  Figure 
1-1)  corresponding  to  the  polarization  p  may  be  expressed  as 

«p<£o>  =  1  -  lRp|2  "  i  J  tVph(  Jsorjk  >  +  Vpv(k0.!sn  sin6d6d4>  (1) 

where  the  integral  is  evaluated  over  the  upper  hemisphere  0  5  6  <  TI/2  , 

0  <  4  <  2r  ,  Here  Rp  is  the  reflection  coefficient  of  the  medium  which,  in 


1-1 
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The  main  problem  to  be  addressed  in  this  section  is  the  computation 
of  the  reflection  coefficient  Rp  and  scattering  coefficients  ^ab  when  the 
random  medium  is  described  by  a  tensor  dielectric  constant  and  correlation 
function  appropriate  to  sea  ice  as  discussed  in  [2].  The  atmospheric  effects 
mentioned  above  may  be  treated  by  standard  radiative  transfer  theory. 

1.2  THE  REFLECTION  COEFFICIENT 

According  to  the  bilocal  approximation  in  strong  fluctuation  theory 
the  reflection  coefficient  Rp  (p  =  h  or  v)  in  (1)  and  (2)  is  determined  by 
the  behavior  of  the  mean  electric  field  which,  in  turn,  is  dependent  on  the 
effective  dielectric  constant.  For  sea  ice  with  no  preferred  azimuthal 
direction,  this  is  a  second  rank  x  msor  of  the  form  [2] 


if  the  incident  wave  is  kQ  =  ko(sin8o,0,-cos6o)  .  Since  we  consider  the 

problem  of  a  medium  located  in  the  region  z  <  0  with  air  above,  (4)  holds 

for  z  <  0  while  for  z  >  0  ,  K  is  simply  the  unit  dyad  I  . 

a  a 

Thus,  given  an  electric  field 

£in  *  E0  p  expti  k0(sin60x  “  cos0oz)]  (5) 

incident  on  the  medium  (see  Figure  1-1)  from  above,  the  mean  field  may  be 
expressed  as 
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Ein  +  RpEoP  exp[  i  kD(  sin60x  +  cos0oz)]  if  z  >  0 


E0  f  ( z  )  exp[ikD  sin0ox] 


if  z  <  0  (6) 


The  second  summand  in  the  first  line  of  (6)  represents  the  field  reflected 
in  the  specular  direction  while  the  second  line  of  (6)  represents  the  depth 
dependence  of  the  field  by  the  vector  f(z)  .  If  a  general  layered 
structure  is  considered,  so  that  K  in  (4)  depends  on  z  ,  an  explicit  form 
for  f  cannot  be  written.  However,  f(z)  and  Rp  are  determined  by  the 
fact  that  Em  is  a  solution  of  Maxwell's  equations  with  certain  boundary 
conditions.  Thus  Em  satisfies 

V  x  V  x  Em  -  k2  K  Em  =  0  (7) 

~  ss  ~ 

The  boundary  conditions  require  the  component  of  Em  tangent  to  the  boundary 
surface  at  z  =  0  be  continuous  together  with  the  tangential  component  of 
V x Em  (which  is  proportional  to  the  magnetic  field).  The  boundary  condi¬ 
tions  lead  to  different  results  for  horizontal  and  vertical  polarization  of 
the  incident  wave.  Thus,  these  cases  are  treated  separately. 

Considering,  first,  the  case  of  horizontal  polarization,  the  sym¬ 
metry  of  the  problem  requires  that  f  in  (6)  be  of  the  form 

f(z)  =  f(z)  h  (8) 

where  f(z)  is  a  scalar  valued  function.  Substitution  into  (7)  shows  that 


f  is  determined  by  the  ordinary  differential  equation 
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k22  “  sin20o]  f  =  0 

dz^ 


(9) 


while  the  boundary  conditions  lead  to  the  equation 


df 


dz 


-  i  kQ  co s0o  f  =  -2i  kD  cos0o  at  z  =  0  (10) 


and  the  relation 


=  f(0)  -  1 


(11) 


It  is  interesting  to  note  that  these  equations  are  the  same  as  those  applic¬ 
able  to  a  scalar  dielectric  constant  [1],  [4]  with  the  single  exception  that 
the  scalar  dielectric  constant  is  replaced  by  the  2,2  component  of  the 
dielectric  tensor  for  the  anisotropic  case. 

The  situation  is  more  complex  for  the  case  of  vertical  polarization. 
In  this  case,  the  symmetry  of  the  problem  requires  that  f(z)  in  (6)  have 
two  components: 


f(z)  = 


A  A 

fx( Z  )x  +  f 2 ( Z  )  Z 


(12) 


Upon  substituting  (12)  and  (6)  in  (7)  it  is  found  that 


fz(z) 


i  sin0, 


d  f. 


43 


k0[ K33  -  sin20o]  dz  K33  -  sin26c 


(13) 


'  -1 
*  1 


v 


.-V. 


'■I 


7771 


m 


1-5 


together  with  the  equation 


Rv  =  1  -  f x( 0 ) /cos0o  (16) 

While  the  derivation  of  equations  (13)  -  (16)  required  no  approxi¬ 
mations  beyond  those  inherent  in  the  bilocal  approximation  of  strong  fluctu¬ 
ation  theory,  they  are  extremely  inconvenient  for  numerical  work  because  the 
function  fx(z)  is  an  oscillatory  function  which  must  be  found  over  dis¬ 
tances  of  many  wavelengths  in  typical  applications  to  snow-covered  sea  ice. 
The  numerical  techniques  discussed  in  [1]  would  be  useful  and  applicable 
except  for  the  possible  non-zero  values  of  and  its  derivative.  For 

this  reason,  these  terms  will  be  ignored  so  that  the  approximate  equations  to 


be  solved  for  vertical  polarization  are 


li 


■«  M  'I 


JTFjmWUT*!.  ■*  ■jm'i  *X  M  fjl  7M  WTT.VT  J  ’X  «_1  «JI  >il 


•y  *y  •}!  *v  u  v  ■)  ■»  'j  'vji'j 
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£_£> 

dz2 


sin26c 


»J3  -  *i"2e„ 


/dK33  /  I  \££x  2  ^11 
\  dz  /  33  /  dx  °  k^3 


(K33  - 


sin^60)  fx 


i  sin0o 


df. 


kofK33  ~  sin260I  dz 


0 

(17) 

(18) 


with  the  boundary  condition 


i 


K33  -  sin260 

ko  1 

K33  cos60 


-2i  k. 


K33  -  sin26c 


^33 


( 19) 


at  z  =  0  .  Here  K33  and  K33  are  slightly  modified  values  of  K33  and 
K33  .  The  discussion  below  treats  the  validity  of  this  approximation  and  the 
fact  that  K31  and  K33  should  be  modified  when  is  ignored. 

With  equations  (8)  -  (11)  for  horizontal  polarization  and  (16)  - 
(19)  for  vertical  polarization,  the  numerical  problems  for  a  multi-layered 
anisotropic  structure  are  exactly  the  same  as  those  treated  for  the  isotropic 
case  in  [1],  Thus,  the  specif ication  of  the  mean  field  and  reflection  coef¬ 
ficients  for  both  horizontal  and  vertical  polarization  may  be  considered  com¬ 
pleted  when  the  components  of  the  dielectric  tensor  are  given.  For  sea  ice, 
these  are  discussed  in  [2]. 

The  nature  of  the  hypothesis  that  K^3  may  be  ignored  for  the  case 
of  vertical  polarization  will  now  be  examined.  To  do  this  notice  that, 
according  to  the  analysis  in  [21,  the  equation  (4)  may  be  written  more 
generally  as 


K  =  AI  +  Bkk  +  Czz  +  D(kz  +  zk) 


(20) 
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where  the  scalar  coefficients  A  ,  B  ,  C  ,  and  D  depend  only  on  k*z  and 
k2  .  For  frazil  or  multiyear  ice  which,  according  to  the  model  in  [2],  are 
characterized  by  an  isotropic  mean  structure,  the  coefficients  C  and  D  are 
exactly  equal  to  zero.  But  then  the  analysis  in  |5)  shows  that  for  trans¬ 
verse  waves  the  coefficient  B  has  no  influence  on  the  effective  dielectric 
constant  so  that  K  may  be  treated  as  equal  to  the  scalar  A  .  Thus,  for 

'■'w 

frazil  and  multiyear  ice,  no  approximation  is  made  at  all.  For  all  types  of 
ice,  no  approximation  is  made  at  60  =  0  because  it  was  shown  in  [2]  that 
*13  =  -  whenever  kx  =  0  .  There  is  an  approximation,  however,  for  ether 
angles  for  anisotropic  ice  forms.  In  order  to  treat  this  case,  consider  a 
uniform  medium  without  gradients  so  that  fx(z)  in  (^2)  has  an  exponential 
z  dependence  exp[-ikzz]  «  exp ( - i kG ( K- s in20o  )  z  )  where  K  is  the  effective 
dielectric  constant  appropriate  to  0O  .  If  it  is  required  that  kz  be 
equal  to  the  exact  propagation  constant  (i.e.,  obtained  as  in  [2]  without 
neglecting  *13'  then  by  substituting  the  exponential  into  ^7)  and  specifying 
tha  t 


K 


11 


K  2  t 


t  21  ) 


it  is  found  that  must  be  given  by 


'33 


s  i  n^Q  K 


K  i  *  s  l  n 


Thus,  a  procedure  has  been  found  for  mo  i  1  f  v  1  n  1  the  die  le>~t  r  1  ~  constants  which 
does  not  change  the  propagation  constant  of  the  wave  when  K  j  5  is  ignored. 
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However,  there  is  a  change  in  the  ratio  f2(z)/fx(z)  which  is  induced  by 
this  approach.  From  (18),  the  approximate  ratio  is 


fz/fx 


sin0o  \/ K-sin20o 

*33  -  sin2e0 


(23) 


while  from  (13)  the  exact  ratio  is 


f  z  /f  x 


»in0o  J K-sin20o  -  K 


13 


K33  -  sin*0c 


(24) 


Numerical  comparisons  for  a  variety  of  anisotropic  ice  parameters  and  angles 
0O  show  that  (23)  yields  values  close  to  (24).  The  largest  differences 
occur  for  high  salinities  and  temperatures  near  0°C  where  differences  of  the 
order  of  10%  in  both  the  real  and  imaginary  parts  of  the  ratios  were  found 
for  large  0O  and  frequencies  near  20  GHz.  At  frequencies  below  10  GHz  and 
above  30  GHz,  the  error  is  significantly  less.  Computations  from  (23)  and 
(24)  are  nearly  identical  for  low  ice  temperatures  and/or  low  salinity.  Thus 
it  is  seen  that  the  approximation  in  the  dielectric  constant  for  vertical 
polarization  of  the  incident  wave  leads  to  acceptably  small  errors  in  the 
ratio  fz/fx  while  preserving  the  magnitude  and  direction  of  the  propagation 
constant  of  the  electromagnetic  wave. 

1.3  SCATTERING  COEFFICIENTS 

The  computation  of  the  scattering  coefficients  for  sea  ice  follows 
the  pattern  discussed  in  [4],  The  essential  difference  is  that  the  scalar 
correlation  function  is  replaced  by  a  fourth  rank  correlation  tensor  describ¬ 
ing  the  anisotropic  inhomogeneities.  This  adds  a  number  of  new  terms  to  the 
scattering  coefficients  Yab  (a,b  •  h  or  v). 
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Specifically,  the  random  electric  field  £r  in  the  bilocal 
approximation  of  strong  fluctuation  theory  obeys  the  equation  (see  [4]  and 
the  required  substitutions  mentioned  in  [6}  to  convert  the  results  of  [4]  to 
strong  fluctuation  theory) 


Vx  V  x  £r  -  k2  K0  £r  =  -  |  Em 

where  K0  is  the  diagonal  quasi-static  dielectric  tensor 


(25) 


Kc  =  diag  (K01,  K01,  K03) 

rv 


(26) 


discussed  in  [2]  for  sea  ice. 


In  (25),  4  is  the  random  second  rank  tensor 


4  =  AL[I  +  S  AL)"1 
~  w  «  »  « 


(27) 


where  S  is  the  coefficient  of  the  delta  function  in  the  Green's  function 
for  the  operator  on  the  left-hand  side  of  (25)  and 


AL  =  k2(K„  -  Kr) 

te  o  ~ 


(28) 


Here,  Kr  is  the  local  dielectric  constant  (i.e.,  corresponding  to  air  or  an 
ice  grain  which  contains  brine  cells). 

Proceeding  as  in  [4],  the  expected  value  jF ( r )£ Jj * ( r ) (where  the 

*  indicates  complex  conjugate)  may  be  formed  in  terms  of  an  integral  over 


the  Green's  function.  After  a  decomposition  which  represents  the  expected 
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value  as  an  angular  spectrum  of  random  waves  E|(k)  propagating  in  the 
direction  k  (see  Figure  1-1),  it  is  found  that 


o  w 

(2tt)-2  E2  k2  cos©  J  dz'  j  dz"  Aifc(z',k)  Ajm(z",k) 

-  CD  -CD 

f  ^  ( 2  '  )  ^n  ^ 2  ^  ^k^mn^ox  -  ^x'  ^y»  2  »z  ) 


(29) 


where  the  summation  convention  is  used  for  repeated  Latin  letter  indices,  the 
tensor  A  is  the  two-dimensional  Fourier  transform  of  the  Green's  function 
exactly  as  in  [4],  f  is  the  vector  (6)  describing  the  mean  electric  field 
in  the  medium,  and  the  fourth  rank  tensor  W'  is  the  Fourier  transform 

o 

^kimn^i'^'2'  >z" ]  =  jdxdy  C^mn(x,y,z'  ,z")  exp[  i  (£x+ny )  1  (30) 

-  CO 

of  the  correlation  tensor 


ck£mn<x'y'z' ’z")  "  (£k£(x 


+  x' 


y ' *  z ' )  4  mn ( x ' •  y ' ’  2 " ) ^ 


(31) 


which  is  an  even  function  of  x  and  y  .  Finally,  the  scattering  coeffi 
cients  are  related  to  (29)  by  [4] 


Vab(*So'’5^  =  (4ttcos6/E^  cos0o)  b^  bj  ^E^(k)  E^*(k)^  (32) 

To  obtain  practical  expressions  for  Yab  •  necessary  to 

specify  the  symmetry  of  the  tensors  A  and  W'  in  more  detail.  In  the 


1-11 
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case  of  the  Green's  function,  the  physical  situation  under  consideration  is 
azimuthally  symmetric.  Just  as  in  [4],  this  implies  that  A  is  of  the  form 


A4>4)  +  cos24>(App-Ad)^)  cos<t>  sin4>  (App-A^xj,)  cos4»Apz' 
cos<$>  sin4>  (App-A^)  A^+sin24>( App_A4>4>^  sin4>Apz 
cos4>  Azp  sin<*>  Azp  A, 


’zz 


(33) 


Further,  since  the  Green's  function  is  a  solution  of  the  vector  wave  equation 
the  relations 


-tanfl  A 


PP 


(34) 


Azz  =  -tan6  Apj 


(35) 


are  found  to  hold  so  that  there  are  only  three  independent  functions  A^  , 
App  ,  and  Apz  which  need  to  be  evaluated.  These  will  be  considered  in  more 
detail  later.  For  the  present,  we  turn  to  the  tensors  and  wijktl  *n 

order  to  be  able  to  reduce  (32)  to  a  manageable  form. 

It  should  be  noted  that  differs  from  C^jj^  defined  in  [21 

for  sea  ice  only  by  the  fact  that  the  complex  conjugate  is  used  in  the  second 
factor  on  the  right-hand  side  of  (31).  Thus,  using  the  corresponding  assump¬ 
tions  regarding  an  exponential  correlation  for  air  bubbles  with  a  correlation 
length  la  and  volume  fraction  va^r  and  grains  of  ice  with  correlation 
length  £g  and  volume  fraction  vg  ,  one  can  immediately  write 
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ci jk£ ^ r2_rl ^  =  t£airij  lairk£  +  £gij  £gk£  ”  ^airij  £gk0  “  £gij  4airk£ J 


vair  ^ 1  ~  vair>  exPl "  1E2"£iI 


+  Bi  j)t£  ( 0 )  exp[- |r2-r1j/£g] 


(36) 


where 


Bijkf^°*  “  1 4gi j  4gki  ~  £gij  £gk£^  vgrain 


(37) 


The  notation  of  [2]  is  used.  The  most  important  immediate  fact  concerning 
(36)  is  that  and,  hence,  its  Fourier  transform  have  the  same 

symmetry  as  (and  j k f ( 0 ) )  in  12]  with  the  exception  that  the  sym¬ 

metry  under  the  interchange  (ij)— *(k£)  is  lost.  That  is,  they  are  invariant 
under  orthogonal  transformations  leaving  the  z-axis  fixed  and  are  symmetric 
in  the  first  two  indices  and  last  two  indices.  Hence,  it  is  possible  to 
write 


wi  jkf 


'1  ^i  j  4  a2tJ,ik  t>jg  4  6i£  &jkl 


4  bl  dij  6 k 3  l‘i 3  4  b216ik  £>j3  t>i’3  4  bi£  ^ j 3  t>k3 


4  6jk  6i3  btl  4  bj{  bi 3  6k3 1  4  b3  6kf  6i3  6j3 


+  C  di3  6j3  6k3  6t’3 


(38) 


where  the  coefficients  are  scalars.  Since  C'  is  an  even  function  of 
x2~xl  and  y2-yi  •  (3°'  and  (31)  show  that 
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wk£ij  =  wijk£  * 

Thus,  using  (38),  it  is  found  that  »  b2  an<^  c  are  real  while 
b|*  «  bj  .  Furthermore,  choosing  particular  values  for  the  indices  shows 
that  the  coefficients  in  (38)  are  related  to  W'  by 


W1 122  * 

1 

al 

w1212  “ 

1 

a2 

W1 1 33  * 

• 

al 

♦  bj 

w1313  ' 

f 

a2 

+  b2 

w  3  3  3  3  * 

1 

al 

+  2a 

aj  +  2aj  +  ^1  +  ^2  +  4b2  +  c 


which,  with  the  exception  of  WH33  *  are  real.  This  set  of  equations  is 
easily  inverted  to  yield 


a  i  =  W 


1122 


w1133  "  w1122 


a2  *  w1212 

b2  *  W1 31 3  ”  W1 2 1 2 

c'  “  w3333  ”  2Re  W1 1 3 3  +  w1122  "  4W1313  +  2W1212 
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of  the  Green's  function  and  correlation  tensor  are  of  the  same  nature  as  in 
11] .  It  will  be  shown  below  that  this  is  indeed  the  case  so  that  computer 
codes  based  on  the  same  numerical  scheme  are  applicable.  Of  course,  they 
will  be  considerably  longer  because  of  the  new  terms. 

1.4  GREEN'S  FUNCTIONS 

The  Green's  function  F(r,r')  to  be  used  in  computing  the  scatter- 

%  ~ 

ing  coefficients  is  defined  as  the  solution  of  the  equation 

V  x  V  x  T-  k2  K0  T  =  I  6  ( r  -  r')  (46) 

«  °  «  a 

where  is  the  quasi-static  dielectric  tensor  (26).  If  the  Fourier 

a;0 

transform 


CD  CD 

Tij  =  (2tt)-2  j  d  4  j  dr)  aj_j(z,z' exp ji[ |(x-x' )+q(y-y ' ) ]j 

-00  -  CD 

(47) 

is  introduced,  then  the  tensor  A  in  (33)  is  defined  in  terms  of  the  value 
at  z  =  0  of  a^j  for  the  special  values  of  the  transform  variables 

£  =  kQ  sinQ  cos4i 

r)  =  kQ  sin0  sin6  (48) 

as  in  [4J.  Explicit  differential  equations  for  a^j  are  found  by  substitut¬ 
ing  (47)  into  (46).  Since  the  symmetry  of  a^j  is  exactly  the  same  as  that 
of  A^j  ,  the  three  independent  functions  describing  A  are  readily  found  to 
be  defined  by  the  following.  For  : 
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A^z*  ,k)  -  a^O) 

where 


d2a<t,(t> 

dz2 


+  ko(K01 


sin2 0)  a^ 


-6(2  -  Z'  ) 


and 


da  <£,<+, 
dz 


i  kocos0  a^cp 


0  at  z  =  0 


For  App  : 


where 


App(z',k)  =  app( 0 ) 


d2at 


sin2  6  dK03  dap0  +  ^2  K01 

6z*  Kg3 ( k02  ~  sin20)dz  dz  +  °  K03 


(k03 


and 


Kq3  -  sin20 
k03 


6  ( z  -  z  '  ) 


dag 

dz 


^“pp  k03  _  sinz0 

-  i  k„  — - - T—  aon  =  0  at  z  =  0 


°  Kq3  cos0  PP 


For  Apz  : 


where 


ApZ ( z ' , k )  =  apz ( 0 ) 


d2c 


sin 


20 


dK03  dapz  ,  2  K01 

~  +  k„  “  (Kq3 


dz2  K03(K03  -  sin20 )  dz  dz  °  k03 


i  sin0 
ko  k03 


d  6(  z  -  z '  ) 
dz 


dK 


03 


Kq3  -  sin 


20  dz 


sin20)app 


sin20)apz 


?.(  z  -  z  '  ) 
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and 


dapz  i  kQ( Kq3  -  sin28) 

dz  Kq3  cosQ  aP2 


i  sin0 
ko  k03 


6  (  z  '  ) 


at  z 


0 


(57) 


In  all  cases  the  boundary  condition  as  z  —  - ®  is  that  A  be  an  outgoing 
wave . 

While  slightly  more  complicated  than  the  corresponding  equations 
for  the  isotropic  case  where  the  dielectric  constant  is  a  scalar,  these 
equations  are  still  of  the  same  mathematical  form  as  treated  in  [1).  Thus, 
equivalent  Riceati  equations  may  be  introduced  and  solved  numerically  as 
discussed  in  [11- 
1-5  CORRELATION  FUNCTION 

The  only  remaining  problem  to  be  discussed  in  connection  *.ith  the 
computation  of  the  scattering  coefficients  is  the  treatment  of  the  Fourier 
transform  of  the  correlation  tensor  (31).  The  first  step  is  the  computation 
of  the  Fourier  transforms  of  the  two  exponential  functions  in  (36).  Lettina 

X  X 

( 4 , p, z , z '  )  =  j  dx  |  dy  exp[ - |r-r ' | /£)  exp[ -i ( | ( x-x ' )  +  q(y-y ' )  (58) 

—  cc  - 

when  £  may  be  either  of  the  correlation  lengths  £a  or  £g  in  (36),  it  is 
seen  that  W^fmn ( kox-kx , ky , z ’ , z " )  in  (29)  has  a  kG,k,  and  spatial  depend¬ 
ence  defined  by  the  functions  ( kox-kx , ky , z , z ' )  ( £  *  £a  or  £g  ).  The 

integral  in  (58)  may  be  evaluated  explicitly  to  yield 

Vp  =  2m£2  u”3[ 1  +  u|z  -  z'|/£]  exp[-u|z  -  z'|/£]  (59) 
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where 

u  =  [1  +  £2(£2  +  q  2  )  ]  1/2  (60) 

Thus,  upon  using  (36)  and  (37),  the  components  of  W  occurring  in  (42)  - 
(45)  nay  be  expressed  as 

w1122  =  Bil22(0)  vi?g  +  1 1  air  "  ^  cgl  1 1 2  vair(1  ~  vair>  v£a  (61) 

W1 2 1 2  =  b1212(0’  vt g  (62) 

W1 1 3 3  =  B1133(0)  v£g  +  (4airll  “  4gll)(4air33  ~  4g335  ^ir'1  "  vair 

(63) 

w1313  =  B1313f0)  vi?g  (64 

w3333  =  b3333^0^  v£g  +  I^air33  “  4g33l2  vair^  _  vair'  v£a  (65) 


The  remainder  of  the  problem  is  identical  to  the  evaluation  of  the 
tensor  in  [2].  Approximations  have  to  be  made  when  the  probability 

of  the  orientation  of  the  ice  grains  has  a  preferred  direction  just  as  in  [2] 
For  frazil  and  old  ice,  where  all  solid  angles  are  equally  probable,  a  more 
exact  computation  is  possible.  Since  the  details  of  the  computation  are 
lengthy,  but  follow  that  in  [2]  (differing  only  in  the  occurrence  of  terms 


containing  complex  conjugates),  only  the  results  will  be  given.  Both  4  air 
and  the  components  of  the  tensor  4„  are  the  same  as  in  [2]  so  that  only  the 


tensor  B'  needs  to  be  described.  Using  the  notation  of  [2],  for  all  cases 
except  for  frazil  and  old  ice. 
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,  (3  S3j32cos6Bsin2esl2  /- 

b1122(0)  '  |  4  axY  I  "  Re^gll 


S3P2cos26Bsin26B 


1  psin26B( P3+S3 ft] I2  j 

i  °i^  1  j 

B'  ..(0)  =  l|psin26B[o  3+s3{3]  2  „ 

1212  g  - 5—^  g 

,  1  pcos6Bsin6g  2 

b1313  f ®  )  =  2  V  vgrain 

Bil33(0)  =  b3333(0)  =  0 


'gram 


gram 


where  0B  is  the  mean  angle  between  the  axes  of  the  brine  cells  and  the 
vertical.  In  the  case  of  frazil  and  old  ice. 


b1212(0) 

=  b{313(0) 

=  J_ 

15 

P 

v  v 

b1122(0) 

=  83333(0) 

■  i 1 

■  Bi j ji 

",3333(0) 

=  Bil22(0) 

+  2B3232 ( 0 ) 

where  the  doubly  contracted  term  in  (71)  is  given  by 


Bijji(0) 


;Tr(P  PM  -  3  Tr  P  ,  vgrain 


Here,  the  second  rank  tensor  P  is  the  diagonal  tensor 
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d 1  aa 


ko( Kof-*l 1 
1  *  k2sf<Kof-Kj) 


ko 1  *of  ~  3  '  \ 

1  ♦  k2Sf'Kof-K3  '/ 


kof*of-*l  1 
1  ♦  k2S{iKQf-Kx) 


f  -’4 


and  Tr  mean  trace. 
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Section  2 

THE  BRIGHTNESS  TEMPERATURE  OP  SEA  ICE 


The  theory  developed  in  Section  1  may  be  applied  to  the  computation 
of  the  brightness  temperature  of  various  types  of  sea  ice.  It  will  be  seen 
that  it  leads  to  predictions  of  microwave  brightness  temperature  signatures 
that  do  not  necessarily  fall  neatly  into  the  categories  which  have  been 
developed  by  the  World  Meteorological  Organization  (WMO)  for  the  visual 
identification  of  ice  types.  A  difficulty  of  this  nature  has  also  been  noted 
in  [1]  where  a  comparison  of  the  ice  types  distinguishable  by  the  radiometric 
observations  discussed  in  [1]  is  shown  not  to  discriminate  between  young  ice 
and  first-year  ice.  A  different  classification  was  proposed  in  [1]  (see 
Table  2-1).  This  classification  still  has  an  element  of  dependence  on  ice 
thickness.  However,  the  computations  based  on  the  theory  of  this  report 
indicate  that  the  deciding  factors  determining  the  microwave  signature  are, 
except  for  the  very  thinnest  of  ice  covers  (of  the  order  of  1  cm  thickness 
where  the  underlying  sea  water  may  influence  the  signature),  thickness 
independent.  The  critical  factors  are  the  ice  salinity  and  the  presence  or 
absence  of  snow  cover  for  new,  young,  and  first-year  ice.  The  conventional 
WMO  categories,  however,  will  be  retained  here  in  order  to  facilitate  compari 
sons  with  published  results.  New,  young,  first-year  and  multiyear  sea  ice 
are  discussed  below. 

2.1  NEW  ICE 

Although  reports  on  the  microwave  brightness  temperature  of  sea  ice 
have  appeared  in  the  literature  since  1972  [2]  and  distinguished  between  "old 
and  "new"  ice,  a  more  systematic  classification  of  the  "new"  ice  category 


Report  8362 


occurred  only  in  later  studies.  Several  recent  observations  of  new  ice  have 
been  reported  in  the  literature.  These  include  those  of  the  NORSEX  Group  [3] 
Hollinger,  et  al  [1],  and  Troy,  et  al  [4].  Before  examining  the  results  of 
these  particular  experiments,  it  will  be  useful  to  make  some  observations  on 
the  physical  characteristics  of  sea  ice  from  its  formation  to  the  beginning 
of  the  melting  season. 

A  very  useful  description  of  these  characteristics  is  given  by 
Nakawo  and  Sinha  [5].  They  observed  ice  at  Pond  Inlet,  Canada  from  freeze-up 
through  May  with  a  careful  documentation  of  salinity  profiles,  temperature, 
and  snow  cover.  An  important  observation  from  the  point  of  view  of  microwave 
radiometric  studies  was  that  the  initial  salinity  of  the  ice  after  formation 
can  be  extremely  high.  For  example,  a  value  of  25°/oo  was  measured  in  one 
case  and  the  authors  commented  that  even  this  may  not  be  the  upper  bound 
because  of  possible  brine  drainage  during  the  removal  of  the  core  from  the 
ice  sheet.  This  initial  salinity  drops  very  rapidly  to  a  more  stable  value. 
In  the  example  cited,  the  salinity  had  dropped  to  8.5°/oo  within  a  week. 
Afterwards,  it  remained  relatively  constant  with  a  tendency  to  decrease 
slowly  during  the  rest  of  the  season.  The  air  temperature  was  generally 
below  -20°C  during  these  measurements.  For  warmer  days,  one  would  expect  the 
high  initial  salinity  to  drop  even  faster  because  of  the  enhanced  brine 
drainage.  The  growth  rate  of  the  ice  depends,  of  course,  on  the  air  tempera¬ 
ture.  Typical  rates  ranged  from  about  0.7  to  1.7  cm/day  for  the  temperatures 
prevailing  at  Pond  Inlet  (-20  to  -40°C  except  during  a  warmer  period  in  the 
first  few  weeks  ).  The  snow  cover,  which  will  be  of  interest  below,  was  also 
monitored.  Within  a  few  weeks  of  the  freeze-up  the  snow  cover  was  already 
about  8  cm  thick  and  throughout  the  remainder  of  the  season  it  fluctuated 
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around  11.4  14.2  cm.  The  snow  density  was  0.35  10.04  g/cm3  during  most  of 
the  season  in  agreement  with  earlier  reported  results.  However,  it  was  men¬ 
tioned  that  the  density  may  have  been  lower  than  this  during  the  early 
growth  period. 

Turning  to  the  microwave  observations,  the  NORSEX  results  [3]  will 
be  considered  first.  The  observations  were  made  between  the  end  of  September 
and  early  October  1979.  The  ice  temperature  was  of  the  order  of  -1°C.  Meas¬ 
urements  at  frequencies  of  4.9,  10.4,  21,  36  and  94  GHz  with  both  vertical 
and  horizontal  polarization  were  made  at  a  nominal  angle  of  incidence 
0  =  50°  (this  varied  by  about  13°  due  to  heeling  of  the  shipboard  platform). 
Figure  2-1  is  a  plot  of  the  deduced  emissivity  (i.e.,  after  removing  atmos¬ 
pheric  effects  and  dividing  by  assumed  ice  temperature  of  272  K)  of  the  calm 
sea  near  the  ice  together  with  newly  frozen  pancake  ice  (case  C  in  Figure  4 
of  [3])  for  both  horizontal  and  vertical  polarization.  The  computed  sea 
water  emissivity  is  also  shown  and  is  mainly  of  interest  as  a  check  on  the 
atmospheric  correction  procedure  and  the  calibration  of  the  instruments  since 
the  radiometric  behavior  of  sea  water  is  relatively  well  understood.  Com¬ 
puted  values,  assuming  an  ice  density  p  «  0.9  g/cm3,  a  salinity  S  * 
12.3°/oo,  frazil  (i.e.,  random  ice  grain  orientations  in  the  theory  of 
Section  1)  with  a  mean  air  bubble  diameter  da  -  0.5  mm  and  a  grain  diameter 
dg  -  4  mm,  and  no  snow  cover,  are  also  shown  for  an  angle  of  incidence  6  « 

50“ .  The  agreement  is  quite  good.  Although  the  salinity  of  the  ice  was  not 
stated  in  [3],  the  assumed  value  appears  to  be  reasonable  on  the  basis  of  the 
warm  temperatures  and  the  discussion  above.  Computations  (not  shown)  were 
also  performed  for  columnar  ice  with  a  mean  brine  cell  angle  from  the  vertical 
of  24°.  In  this  case,  a  salinity  of  13°/oo  (which  is  still  quite  reasonable) 
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had  to  be  assumed  in  order  to  obtain  a  fit  to  the  experimental  data  which  is 
roughly  as  good  as  that  shown  for  frazil.  For  columnar  ice,  S  «  12.3°/oo 
yields  a  distinctly  poorer  fit  to  the  data.  In  view  of  the  fact  that  the 
uppermost  layer  of  ice  is  generally  frazil  or  randomly  oriented  columnar  ice 
[II,  [61,  the  initial  assumption  is  probably  more  realistic  although  the 
microwave  evidence  does  not  indicate  a  distinct  preference  for  one  or  the 
other  of  these  models.  The  underlying  sea  water  has  almost  no  effect  on  the 
computed  emissivity  for  an  ice  thickness  of  1  cm  because  of  the  high  attenua¬ 
tion  through  the  ice.  Thus,  for  all  practical  purposes,  the  ice  thickness 
has  no  direct  bearing  on  the  emissivity.  A  characteristic  feature  of  these 
results  is  the  monotonic  increase  of  the  emissivity  with  increasing  frequency 
This  can  be  traced  directly  to  the  fact  that  at  this  high  salinity  and  high 
temperature,  both  of  which  imply  a  large  brine  volume,  the  dielectric  con¬ 
stant  of  the  sea  ice  is  quite  large  at  the  lower  frequencies  but  gradually 
decreases  as  the  frequency  increases. 

For  the  warm  ice  temperatures  under  the  NORSEX  conditions  one  would 
expect  rapid  surface  desalinization  and,  perhaps,  a  lowering  of  the  ice  sur¬ 
face  density  due  to  brine  drainage  as  the  ice  develops.  Thus,  it  is  not 
surprising  that  for  the  densely  packed  3  cm  thick  pancake  ice  (case  D  in 
Figure  4  of  [3])  the  emissivity  is  considerably  higher  than  for  the  thin 
pancake  ice.  Figure  2-2  is  a  plot  of  this  data  together  with  computations 
which  assumed  an  ice  density  p  «  C.8  g/cm3,  salinity  S  -  9.3°/oo,  mean  air 
bubble  diameter  dfl  *  0.5  mm,  and  mean  grain  size  dg  «  4  mm.  Again,  no 
preferred  orientation  of  the  ice  qrains  is  assumed  although,  with  an  adjust¬ 
ment  of  the  salinity,  similar  results  are  obtained  for  columnar  ice  with  a 
mean  brine  cell  angle  of  24°.  It  is  seen  that  while  the  computations  for 
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vertical  polarization  agree  well  with  the  measurements,  they  are  too  low  for 
horizontal  polarization.  There  does  not  seem  to  be  a  ready  explanation  for 
this  discrepancy.  As  in  the  case  of  thin  pancake  ice,  both  theory  and  experi¬ 
ment  show  the  same  characteristic  increase  of  emissivity  with  increasing 
f  requency . 

Another  set  of  data  on  snow-free  new  ice  is  found  in  [1]  where 
nadir-viewing  (6  *  0°)  measurements  were  made  at  frequencies  of  19,  22,  31, 
and  90  GHz  during  the  freeze-up  period  in  October  1981.  The  ice  temperature 
was  approximately  270  K  during  the  experiments.  Figure  2-3  is  a  comparison 
of  the  experimental  data  (together  with  an  indication  of  the  range  of 
observed  values)  with  computations  which  assume  randomly  oriented  ice  grains 
in  ice  with  p  =  0.85  g/cm3,  da  =  0.5  mm,  dQ  *=  4  mm  and  a  number  of  different 
salinities.  These  parameters  were  not  measured  in  the  aircraft  measurements 
program.  However,  increasing  the  density  to  0.9  g/cm3  reduces  the  computed  emis- 
sivities  slightly.  It  is  seen  that  the  computations  generally  fall  within  the 
measured  range  but  that  the  higher  frequency  data  fall  on  the  lower  salinity 
curves.  This  is  quite  plausible  since  the  smaller  penetration  depth  of  the 
higher  frequencies  would  emphasize  conditions  closest  to  the  surface  where 
some  brine  has  probably  drained  out.  No  attempt  was  made  to  model  this  situ¬ 
ation  because  there  is  no  data  available  on  salinity  gradients  within  a 
centimeter  of  the  surface.  It  should  be  noticed  that  an  increase  of  emis¬ 
sivity  with  increasino  frequency  occurs  in  this  data  just  as  in  the  NORSEX 
experiments  although  the  magnitude  of  the  change  is  smaller.  This  is  due  to 
the  lower  brine  volumes  resultinq  from  the  lower  temperature. 

Data  on  new  ice  at  much  lower  temperatures  is  given  in  [4)  where 
nadir  viewing  (6*0°)  measurements  at  14,  19,  31,  and  90  GHz  were  made  durino 

2-5 


Report  8362 


April  1977.  The  ice  temperature  ranged  between  261  and  264  K.  An  assumed  ice 
temperature  of  263  K  was  used  in  the  computations  shown  in  Figure  2-4.  It  is 
seen  that,  except  for  90  GHz,  the  data  tend  to  follow  a  curve  with  high  density 
and  very  high  salinity.  This  is  entirely  reasonable  considering  the  time  of 
year.  The  ice  will  be  classified  as  new  only  for  a  short  period  because  of 
the  rapid  growth  to  be  expected  so  that  possibly  very  high  salinities  could 
arise  as  pointed  out  above.  The  exception  at  90  GHz  could  be  accounted  for 
by  a  desalinization  of  a  very  thin  surface  layer.  Again,  there  is  a  charac¬ 
teristic  increase  in  emissivity  with  increasing  frequency  which  is  due  to 
the  decrease  of  the  dielectric  constant  with  increasing  frequency. 

Before  leaving  the  subject  of  new  ice,  some  observations  will  be 
made  on  a  possible  model  for  new  and  young  ice  which  has  been  mentioned 
several  times  in  the  literature  (see  e.g.,  [1],  [3],  [4]).  It  is  suggested 
that  these  types  of  ice  have  a  thin  film  of  water  on  the  surface.  Model 
computations  have  been  made  for  new  ice  covered  with  a  film  of  brine  with 
salinity  high  enough  to  prevent  freezing  and  with  thickness  ranging  from 
0.007  to  0.1  mm.  The  results  are  shown  in  Figure  2-5  at  nadir  for  p  *  0.9 
g/cm^,  s  =  10°/oo  and  Tg  =  270  K  and  in  Figure  2-6  for  p  *  0.9  g/cm^ ,  S  ■ 
10.3°/00,  and  Tg  =  263  K.  These  should  be  compared  with  Figures  2-3  and  2-4 
respectively.  Considering  Figure  2-5  first,  it  is  seen  that  a  film  of  brine 
with  a  thickness  of  the  order  of  4  to  20  micrometers  would  be  required  to 
explain  the  observations  below  31  GHz  but,  for  the  assumed  ice  salinity  of 
10°/oo  cannot  explain  the  measurements  at  90  GHz.  However,  just  as  for  the 
computations  of  Figure  2-3,  if  it  is  assumed  that  a  lower  ice  salinity  exists 
near  the  surface  layer  which  affects  90  GHz,  the  computed  emissivity  at  90 
GHz  can  be  raised  slightly.  For  S  =  6°/oo,  the  emissivity  at  90  GHz 
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becomes  0.900  for  a  film  thickness  of  4  micrometers  but  only  0.888  for  a  film 
thickness  of  10  micrometers.  This  is  still  significantly  below  the  measured 
value  and  below  the  value  0.912  achieved  at  90  GHz  for  the  film-free  case  in 
Figure  2-3.  It  is  not  known  if  the  required  film  thickness  below  40  GHz  are 
reasonable.  No  measurements  exist,  but  the  required  values  do  seem  to  be 
small.  Similar  comments  may  be  made  concerning  Figure  2-6  where  a  film 
thickness  of  the  order  of  20  to  50  micrometers  is  required  to  explain  the 
data  below  31  GHz  but  causes  difficulty  at  90  GHz.  Upon  reducing  S  to 
6°/oo  not  much  improvement  is  obtained  because  the  low  temperature  already 
implies  a  low  brine  volume.  Increasing  the  salinity  to  20°/oo  reduces  the 
required  film  thickness  to  about  20  micrometers  up  to  31  GHz  but  again  results 
in  an  emissivity  value  which  is  much  too  low  at  90  GHz.  In  view  of  the  fact 
that  high  frequency  data  in  all  cases  would  be  difficult  to  explain  if  a  film 
were  present,  it  is  probable  that  this  model  is  not  correct.  Rather  new  and 
young  ice  will  feel  moist  to  the  touch  simply  because  the  high  salinity  will 
lead  to  high  brine  volumes  in  the  ice.  This  is  a  volume,  not  a  surface 
effect. 

2.2  YOUNG  ICE 

Both  [1]  and  [4]  present  data  on  young  ice  observed  at  nadir  from 
an  aircraft.  The  emissivity  behaves  differently  in  the  two  experiments. 

While  the  emissivity  data  in  [4]  basically  shows  the  monotonic  increasing 
behavior  with  increasing  frequency  discussed  for  new  ice,  the  data  in  [1]  has 
a  somewhat  higher,  almost  frequency  independent  emissivity  that  is  indistin¬ 
guishable  from  that  of  first  year  ice  (to  be  discussed  below).  While  the  ice 
temperatures  differed  in  the  two  cases,  this  does  not  appear  to  be  the  source 
of  the  differences.  Rather,  the  presence  of  a  thin  layer  of  wet  snow  in  the 
ice  observed  in  [1]  seems  to  account  for  the  difference. 
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Figure  2-7  is  a  plot  of  the  relevant  data.  Measurements  from  [1] 


thickness  ranged  from  1  to  3.5  cm.  Details  of  the  volume  fraction  of  water 
in  the  snow  and  its  density  were  not  published  in  [1]  but  it  was  indicated 
that  at  several  sites  the  bottom  layer  was  slush  while  at  other  sites  it  was 
not  fully  saturated.  For  the  computations,  it  was  assumed  that  volume  frac¬ 
tions  of  water  vw  =  2  and  3%  (far  from  being  slush)  existed  in  the  bottom 
layer  and  that  the  snow  density  was  0.35  g/cm^.  Considering  the  time  of  day 
of  the  measurements,  the  top  layer  could  have  had  a  very  small  amount  of  free 
water.  It  was  assumed  that  vw  ranged  from  0.5%  to  0.8%  in  this  layer. 
Computations  were  averaged  for  a  number  of  total  and  bottom  layer  thicknesses 
in  the  stated  range  and  assuming  a  mean  snow  grain  diameter  of  1  mm.  As  for 
the  underlying  ice,  the  grains  were  assumed  to  have  no  preferred  orientation. 
The  other  ice  parameters  were  p  =  0.85  g/cm^  ,  s  =  10.9°/oo,  mean  air  bubble 
diameter  da  =  0  '  mm,  and  mean  grain  size  dg  =  4  mm.  It  is  seen  that  the 
average  of  the  computations  yield  a  generally  high  emissivity  which  deviates 
by  at  most  0.01  from  a  mean  value  of  slightly  more  than  0.95.  There  is 
little  systematic  frequency  dependence  in  the  computations.  Of  course,  by 
changing  the  snow  and  ice  parameters  somewhat,  one  could  obtain  slightly  dif¬ 
ferent  average  values  than  those  shown.  One  should  compare  these  results 
with  Figure  2-3  which  shows  roughly  the  same  ice  conditions  but  without  snow 
cover.  The  effect  of  the  snow  is  to  raise  the  emissivity  considerably  and  to 
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flatten  the  frequency  response  curve.  The  former  effect  of  snow  on  ice  has 
often  been  pointed  out  in  the  literature.  For  the  snow-free  young  ice  dis¬ 
cussed  in  [4],  the  ice  temperatures  ranged  from  256  to  260  K.  A  -alue  of 
258  K  was  assumed  in  the  computations  together  with  the  same  ice  parameters 
as  above  except  that  the  salinity  was  allowed  to  vary.  The  computed  results, 
which  tend  to  reproduce  the  data,  are  considerably  lower  than  for  the  270  K 
data  and  look  very  much  like  the  computations  shown  for  new  ice  in  Figures 
2-3  and  2-4.  The  main  difference  is  a  shift  to  higher  emissivities  which  is 
readily  understandable  on  the  basis  of  the  lower  ice  temperature  and  hence 
lower  brine  volume. 

2.3  PIRST-YEAR  ICE 

There  are  two  data  sets  for  first-vcar  ice  that  will  be  examined 
here.  First,  consider  the  NORSEX  results  |31  where  the  ice  temperature  was  a 
warm  272  K,  but  the  first-year  ice  was  covered  w : t '  snow.  As  was  mentioned 
in  the  discussion  of  new  ice.  both  vertical  1'.  and  *.  r  :  ?  or  t  a  1  1  y  polarized 
brightness  temperatures  were  measure!  for  an  ancle  f  incidence  fi  =  50® .  The 
NORSEX  data  for  first-year  ice  is  shown  in  Fioure  2-c  where,  except  for  the 
4.9  GHz  and  to  a  much  smaller  extent  the  1'.4  GHz  data,  not  much  frequency 
dependence  is  evident  in  the  emissivity.  Model  computations  for  the  two 
polarizations  were  averaged  over  8,  10,  and  12  cm  depths  of  snow  (a  ranae  was 
chosen  since  the  snow  depth  was  not  specified)  and  a  small  ranae  of  snales 
near  50°  to  simulate  the  angular  uncertainty  in  the  NORSEX  data.  This  proce¬ 
dure  also  had  the  advantage  of  averaginq  over  some  angular  structure  which 
arose  in  the  models  due  to  the  plane  parallel  layered  nature  of  the  snow 
model  being  investigated.  Of  course,  in  actual  experiments,  one  does  not 
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expect  to  find  layered  structures  whose  thickness  over  a  radiometric  foot¬ 
print  is  constant  to  within  a  fraction  of  a  centimeter.  As  for  the  snow 
conditions,  one  should  expect  some  moisture  in  the  snow  layer  nearest  the  ice 
(recall  the  ground  truth  measurements  for  warm  conditions  mentioned  for  Mould 
Bay  in  the  previous  section)  both  because  of  the  warm  temperature  and  the 
fact  that  salt  from  the  highly  saline  ice  will  infiltrate  the  nearest  snow 
layer  thus  allowing  free  water  to  form  in  the  snow.  Hence,  computations  were 
performed  where  the  bottom  1  or  2  cm  of  snow  were  wet  with  a  volume  fraction 
of  9  to  10%  free  water.  However,  the  bulk  of  the  snow  ( 7  to  9  cm  thick)  was 
allowed  to  be  either  dry  or  with  a  very  small  amount  of  free  water  (0.5%)  to 
simulate,  perhaps,  a  trace  of  meltinq  due  to  solar  radiation,  in  all  cases, 
the  snow  density  was  assumed  to  be  0.35  g/cm^ .  The  ice  conditions  were 
assumed  to  be  p  =  0.85  g/cm^ ,  s  =  10.9°/oo,  da  =  0.5  mm,  dg  =  4  mm  and  no 
preferred  orientation  for  the  grain  axes  in  the  uppermost  part  of  the  ice. 
Computations  for  the  two  assumed  snow  conditions  are  shown  in  Figure  2-8. 
Considering  the  uncertainties  in  the  assumed  conditions,  it  is  seen  that  the 
computations  are  reasonably  close  to  the  measurements  with  the  case  where  the 
upper  snow  layer  has  a  trace  of  free  water  being  slightly  better  than  the  dry 
upper  layer. 

Much  colder  first-year  ice  (Tg  =  25C  to  257  K)  was  measured  at 
nadir  in  the  aircraft  program  described  in  [4].  The  data  are  shown  in  Figure 
2-9  together  with  the  results  of  computations  with  two  different  snow  cover 
models,  both  at  a  temperature  Tg  =  253  K.  The  lower  curve  represents  average 
values  for  dry  snow  with  a  range  of  thicknesses  up  to  9.5  cm,  a  density  of 
0.35  g/cm^  and  mean  grain  sizes  of  0.7  and  1  mm.  The  upper  curve  represents 
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an  average  for  a  layered  snow  model  with  a  volume  fraction  of  free  water 
vw  =  1  or  2%  in  the  bottom  layer  of  thickness  ranging  from  2  to  2.5  cm  to 
simulate  possible  infiltration  of  salt  into  the  snow  and  an  upper  layer  with 
a  trace  of  free  water  (vw  =  0.5%)  to  represent  possible  effects  of  solar 
radiation.  The  total  snow  thickness  for  this  situation  ranged  from  9  to  12.5 
cm  and,  again,  the  density  was  taken  to  be  0.35  g/cm3.  The  ice  in  both  cases 
was  assumed  to  have  an  upper  layer  of  frazil  with  p  =  0.85  g/cm3,  da  =  0.5  mm 
and  dg  =  4  mm.  For  the  dry  snow  over  ice  cases,  the  ice  salinity  was  10.9 °/oo, 
but  was  reduced  slightly  to  10°/oo  for  the  second  set  of  snow  conditions.  The 
model  with  a  trace  of  water  yields  emissivities  closer  to  the  measurements,  but 
there  is  no  independent  information  to  support  the  assumptions  made.  Consider¬ 
ing  the  expanded  scale  and  lack  of  definitive  knowledge  of  snow  conditions,  the 
agreement  of  the  measurement  and  computations  in  Figure  2-9  is  reasonably  good. 
The  principal  cause  for  possible  concern  is  the  tendency  for  the  computed  emis¬ 
sivities  to  fall  somewhat  faster  with  frequency  than  the  measurements  above 

20  GHz.  A  possible  reason  for  this  is  the  fact  that  the  loss  tangent  of  the 

ice  grains  in  the  snow  at  the  low  temperatures  prevailing  here  is  not  well 
known  at  these  frequencies.  Since  this  has  an  influence  on  the  computed  scat¬ 
tering  from  the  snow,  the  computations  are  somewhat  uncertain.  Increasing 
the  loss  tangent  would  raise  the  computed  emissivity  at  high  frequencies. 

2.4  MULTIYEAR  ICE 

The  emission  characteristics  of  multiyear  ice  have  been  discussed 
since  the  earliest  microwave  observations  of  sea  ice.  Usually,  no  finer  dis¬ 
tinctions  concerning  the  ice  age  is  made  but  occasionally,  as  in  [1),  second- 
year  ice  is  separated  from  older  ice.  In  the  present  discussion,  all  multi¬ 
year  ice  is  treated  as  one  kind  of  ice  not  only  because  of  a  lack  of 


2-11 


Report  8362 


microwave  data  over  a  wide  range  of  conditions  but  also  because  a  clear 
delineation  of  differences  in  the  physical  structure  of  various  age  subdivi¬ 
sions  is  not  available. 

First  consider  relatively  cold  multiyear  ice  such  as  observed  in 
[4]  in  the  temperature  range  247-256  K  at  nadir.  This  may  be  supplemented 
by  the  data  of  Gloersen,  et  al  [7]  for  ice  at  250  K,  Campbell,  et  al  [8]  for 
ice  in  the  257-260  K  range,  and  Gloersen,  et  al  [9]  for  ice  near  260  K.  Figure 
2-10  is  a  plot  of  the  measured  nadir  emissivities  in  the  frequency  ranqe  below 
40  GHz.  No  indication  of  the  variability  was  given  in  [7]  and  [9],  but  it  is 
reasonable  to  expect  the  same  variability  as  indicated  in  [4]  and  [8].  Snow 
was  always  present  on  the  multiyear  ice  but  details  were  not  given  in  any  of 
the  reports.  For  the  computations,  ice  with  two  basic  types  of  snow  cover 
were  investigated.  The  first  was  a  uniform  dry  snow  layer  in  the  14-17  cm  thick 
ness  range  with  a  density  of  0.35  g/cm3  and  mean  grain  diameter  of  1  mm.  In 
some  cases  a  1  to  2  mm  thick  crust  with  a  density  of  0.5  g/cm  3  and  1  mm.  mean 
grain  diameter  was  added.  The  second  type  of  snow  cover  was  a  multilayered 
snow  structure  based  on  the  description  of  snow  cover  given  in  [10].  A 
covering  of  snow  with  two  random  fresh  water  ice  layers  embedded  in  the  snow 
was  assumed.  Of  the  14  computed  cases,  the  top,  middle,  and  bottom  snow 
layers  had  mean  thicknesses  and  standard  deviations  from  the  mean  of  4.31 
40.12,  3.07  ±0.27  and  8.90  ±1.3  cm  while  the  two  ice  layers  had  mean  thick¬ 
nesses  and  standard  deviations  of  7.35  ±0.98  and  4.37  ±0.96  mm  leading  to  a 
mean  snow-ice  cover  thickness  of  17.5  cm  on  the  sea  ice.  The  snow  had  a 
density  of  0.35  g/cm3  and  mean  grain  diameter  of  1  mm  while  the  fresh  water 
ice  had  a  density  of  0.9  g/cm3  with  a  0.1  mm  mean  air  bubble  diameter.  The 
old  sea  ice  in  both  the  random  and  nonrandom  snow  cover  cases  was  assumed  to 
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have  a  density  p  «  0.6  q  near  the  surface,  salinity  S  -  0.01°  'oo,  a 

mean  air  bubble  diameter  da  -  1  mm,  and  a  mean  grain  size  dg  ■  4  mm.  The 
temperature  was  set  equal  to  253  K  for  all  of  the  computations.  Average 
emissivity  values  for  the  two  distinct  snow  cover  models  are  shown  in  Fiqure 
2-10.  The  randor  layerinq  leads  to  a  lower  emissivity  than  the  uniform  snow 
cover  model.  Also,  a  detailed  examination  of  the  cases  in  the  first  snow 
model  (not  shown  separately  in  the  fiqure)  shows  that  a  crust  on  the  snow 
tends  to  lower  the  emissivity  slightly.  Considering  the  fact  that  the  true 
physical  state  of  the  snow  cover  is  unknown,  the  computations  are  seen  to 
reproduce  the  main  features  of  the  measurements  to  a  frequency  of  about 
30  GHz.  Above  30  GHz,  the  computations  begin  to  show  an  increase  with  fre¬ 
quency  and  at  90  GHz  average  emissivity  values  of  0.915  and  0.860  for  the 
uniform  and  random  snow  cover  models  are  found.  These  are  somewhat  hiqher 
than  the  value  0.81  ♦  0.04  reported  in  [4].  Although  the  computations  appear 
to  be  too  high  at  90  GHz,  it  is  interesting  to  note  that  particular  samples 
of  data  discussed  in  [4]  show  that  there  may  be  an  increase  in  the  emissivity 
of  multiyear  ice  between  31  and  90  GHz.  Thus,  at  least  qualitatively  the 
upturn  predicted  by  the  theory  at  high  frequencies  is  confirmed.  However, 
the  results  in  [1],  which  are  considered  next,  indicate  that  the  increase  in 
emissivity  does  not  begin  until  a  much  higher  frequency  is  reached. 

The  old  ice  in  [1],  observed  at  a  temperature  Tg  ~  270  K,  was 
characterized  as  second-year,  multiyear  and  old  shorefast  ice.  The  old  shore- 
fast  ice  had  emissivity  values  lying  between  the  second-year  and  multiyear 
ice  and  will  not  be  considered  further  here.  Figure  2-11  is  a  plot  of  the 


measurements  and  shows  that  the  multiyear  ice  has  a  consistently  lower 
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emissivity  than  the  second-year  ice.  Computations  using  exactly  the  same 
sea  ice  and  snow  conditions  as  discussed  for  Fiqure  2-10  except  that  the 
thermal  temperature  is  raised  to  Tg  «  270  K  are  also  shown  in  Figure  2-11. 
Again,  the  true  state  of  the  snow  cover  is  not  known.  Althouqh  the  frequency 
range  of  the  measured  data  is  not  as  large  as  shown  in  Figure  2-10,  one  can 
see  a  general  agreement  between  the  computations  and  the  data  at  least  to  a 
frequency  of  22  GHz.  The  computations  are  too  high  at  31  GHz  for  the  multi¬ 
year  ice  but  not  for  the  second-year  ice.  At  90  GHz,  the  computed  emissivi- 
ties  are  definitely  hiqher  (0.927  for  the  uniform  snow  cover)  than  the  meas¬ 
ured  values  for  either  the  second  or  multiyear  ice  types  (0.84  and  0.67 
respective ly ) , 

One  may  conclude  from  these  examples  that  the  theory  is  able  to 
account  for,  in  a  quantitative  fashion,  the  general  decrease  in  the  emis- 
sivity  of  multiyear  ice  to  a  frequency  somewhere  in  the  20-30  GHz  range. 

The  theory  overestimates  the  emissivity  of  multiyear  ice  at  higher  fre¬ 
quencies.  There  are  two  probable  reasons  for  this  discrepancy.  First,  the 
effect  of  the  snow  cover  is  not  being  computed  correctly  because  the  lack 
of  knowledge  of  the  true  loss  tangent  of  dry  snow.  This  is  more  important 
than  in  the  case  of  first-year  ice  because  of  the  lower  emissivity  of  the 
underlying  multiyear  ice  compared  to  first-year  ice.  Second,  and  probably 
more  important,  is  the  fact  that  the  theory  requires  an  expansion  in  a  small 
parameter.  Truncating  the  expansion  after  the  leading  scattering  terms,  as 
in  the  bilocal  approximation,  underestimates  the  effect  of  scattering  and 
hence  overestimates  the  emissivity  of  the  underlying  sea  ice  (and,  perhaps, 
snow  cover).  The  agreement  between  the  computations  and  measurements  at 
lower  frequencies  is  a  result  of  the  fact  that  scattering  is  less  important 
at  lower  frequencies  so  that  the  truncation  is  valid. 
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TABLE  2-1 .  SEA  ICE  CLASSIFICATION 


Selected  WMO  Terminology 

Microwave 

General 

Specific 

Ice  Thickness 

Classification 

New  Ice  | 

j  f razil , grease , slush , 

<5  -  10  cm | 

new  ice 

|shuga,ice  rindfnilas 

1 

Young  Ice  | 

gray 

10-15  cm  , 

1 

I 

.  gray-white 

15-30  cm  1 

first-year  ice 

. 

( thin 

/ 

30-70  cm  [ 

First-Year 

l  .. 

.medium 

70-120  cm  I 

Ice  j 

i 

l  thick 

>120  cm  ' 

Old  Ice  | 

j  second-year  ice 

up  to  2 . 5  m  or  more 

second-year  ice 

1 

(multiyear  ice 

up  to  3  m  or  more 

multiyear  ice 
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Figure  2-5.  Emissivity  of  Hew  Ice  Covered  with 
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Figure  2-11.  Eeissivity  of  Multiyear  Ice  aa  a  Function  of 
Frequency  for  6-0*  (T^  “  270  K) 
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Section  3 

CONCLUSIONS  AND  RECOMMENDATIONS 

A  theory  of  the  microwave  brightness  temperature  of  sea  ice  based 
on  the  bilocal  approximation  in  strong  fluctuation  theory  has  been  developed 
and  shown  to  be  able  to  account  for  reported  measurements  on  new,  younq,  and 
first-year  sea  ice  over  a  variety  of  temperature  conditions  and  frequencies 
up  to  140  GHz.  Experimental  data  on  multiyear  ice  are  explainable  at  low 
frequencies,  but  the  theory  begins  to  become  less  quantitatively  correct  in 
the  20-30  GHz  range.  Calculations  for  multiyear  ice  at  higher  frequencies 
do  not  accurately  reproduce  measurements. 

The  theory  shows  that  snow  cover  on  sea  ice  has  a  significant 
effect  on  its  emission  characteristics  and  is  responsible  for  the  differences 
in  the  observed  spectral  signature  of  new  and  young  ice  without  snow  cover 
compared  to  young  ice  with  snow  cover  and  first-year  ice  (which  almost  always 
is  covered  with  snow).  The  effect  of  snow  is  not  only  to  increase  the  emis- 
sivity  of  these  types  of  sea  ice  but  also  to  flatten  the  spectral  curve  so 
that  the  emissivity  becomes  almost  independent  of  requency.  For  these  rea¬ 
sons,  young  ice  (which  may  or  may  not  have  snow  cover)  can  appear  to  be  radi- 
ometrically  either  like  new  or  first-year  ice.  This  has  been  observed  in  the 
literature.  Also,  for  these  ice  types,  where  the  salinity  is  high,  it  appears 
that  salt  from  the  ice  can  infiltrate  the  bottom-most  layer  of  snow  and  cause 
the  appearance  of  some  free  water  in  the  snow  even  though  the  temperature  may 
be  quite  low. 

For  snow-free  ice,  a  model  with  a  film  of  water  on  the  ice  surface 
was  investigated.  This  model  has  been  suggested  in  the  literature  to  explain 
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certain  features  of  the  microwave  signature  (which  can  also  be  explained  on 
the  basis  of  a  higher  bulk  salinity).  It  was  found  that  the  required  film 
thickness  is  very  small  -  no  more  than  20  to  30  micrometers.  This  is,  per¬ 
haps,  too  small  to  be  reasonable  so  that  such  a  model  may  not  be  plausible. 

In  addition,  it  does  not  reproduce  hiqh  frequency  (90  GHz)  data  as  well  as 
models  without  the  film.  Because  the  explanation  of  published  data  does  not 
actually  require  such  a  hypothesis,  it  was  rejected  in  this  study. 

The  principal  difficulty  with  the  theory  appears  to  be  in  account¬ 
ing  for  the  high  frequency  (above  30  GHz)  decrease  in  the  emissivity  of  multi 
year  sea  ice.  A  possible  cause  for  this  problem  may  lie  in  the  fact  that  the 
loss  tangent  of  pure  ice  (as  it  exists  in  the  ice  grains  of  the  snow  cover) 
is  incorrectly  specified  and  thus  leads  to  a  poor  estimate  of  scattering 
effects  in  the  snow  and  its  consequent  influence  on  the  ice  emissivity.  The 
lack  of  knowledge  concerning  the  dielectric  constant  of  pure  ice  has  been 
pointed  out  in  previous  studies.  Another,  possibly  more  significant,  reason 
for  the  failure  to  predict  the  small  emissivity  of  multiyear  sea  ice  at  high 
frequencies  is  the  fact  that  the  bilocal  approximation  in  stronq  fluctuation 
theory  ignores  second  and  higher  order  terms  in  a  perturbation  expansion  of 
the  random  field  scattered  by  the  inhomoqeneities  of  the  medium.  The  reduced 
scattering  implies  a  higher  emissivity.  At  lower  frequencies,  where  scatter¬ 
ing  may  be  adequately  treated  by  keeping  only  the  leading  term,  the  theory 
was  able  to  reproduce  the  experimental  data  quite  well.  Thus,  the  develop¬ 
ment  of  a  method  for  including  higher  order  scattering  corrections  should 
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On  the  experimental  side,  this  study  has  shown  that  ice  conditions 
such  as  density  and  sa  1  ini ty  qradien t s  within  the  first  few  centimeters  of  the 
surface  (or  ice-snow  boundary)  may  be  quite  important.  Yet,  almost  no  experi¬ 
mental  data  is  available  on  this  subject.  Much  more  attention  should  be  paid 
to  this  area  in  studies  related  to  microwave  radiometry.  Also,  the  question 
of  whether  future  models  should  consider  possible  water  films  on  the  ice 
depends  on  measurina  a  film  of  water  only  a  few  tens  of  micrometers  thick. 
Tentatively ,  this  node  1  has  been  reiected,  but  a  definitive  decision  requires 
more  data.  Finally,  an  extremely  important  area  from  the  point  of  view  of 
radiometric  studies  is  a  full  description  of  any  snow  cover  on  the  ice. 
Descriptions  of  this  nature  have  almost  never  appeared  but,  as  this  study  has 
shown,  the  snow  cover  is  very  significant.  Whether  it  is  dry  or  wet,  -jr.  i  for- 
o-  multi  layered  with  embedded  ice  layers,  thick  or  thin  can  have  an  important 
role  in  the  interpretation,  com;  utat  l  or. ,  and,  ultimate  achievement  of  a  full 
understanding  of  sea  i -e  emissivity.  Of  course,  to  obtain  such  da*a  renui res 
ground  truth  data  tha*  is  often  not  available  when  aircraft  or  satellite  over¬ 
flight  of  i cf  are  made.  However,  as  at  least  one  example  in  this  st  idv  has 
shown,  such  data  car  be  very  valual  le  so  tha*  occasional  supporting  'ir-’in  ! 
parties  should  be  fielded  to  surport  these  remote  experiments. 


